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We studied the properties of the large-scale circulation (LSC) in turbulent Rayleigh-Benard (RB) 
convection by using results from direct numerical simulations in which we placed a large number 
of numerical probes close to the sidewall. The LSC orientation is determined by either a cosine or 
a polynomial fit to the azimuthal temperature or azimuthal vertical velocity profile measured with 
the probes. We study the LSC in F = D/L = 1/2 and F = 1 samples, where D is the diameter 
and L the height. For Pr — 6.4 in an aspect ratio T = 1 sample at Ra = 1 X 10 8 and 5 x 10 8 the 
obtained LSC orientation is the same, irrespective of whether the data of only 8 or all 64 probes per 
horizontal plane are considered. In a T = 1/2 sample with Pr — 0.7 at Ra — 1 x 10 8 the influence 
of plumes on the azimuthal temperature and azimuthal vertical velocity profiles is stronger. Due 
to passing plumes and/or the corner flow the apparent LSC orientation obtained using a cosine fit 
can result in a misinterpretation of the character of the large-scale flow. We introduce the relative 
LSC strength, which we define as the ratio between the energy in the first Fourier mode and the 
energy in all modes that can be determined from the azimuthal temperature and azimuthal vertical 
velocity profiles, to further quantify the large-scale flow. For Ra = 1 x 10 s we find that this relative 
LSC strength is significantly lower in a V = 1/2 sample than in a T = 1 sample, reflecting that the 
LSC is much more pronounced in a V — 1 sample than in a V = 1/2 sample. The determination 
of the relative LSC strength can be applied directly to available experimental data to study high 
Rayleigh number thermal convection and rotating RB convection. 



I. INTRODUCTION 

Rayleigh-Benard (RB) convection is the motion of 
a fluid contained between two parallel plates which is 
heated from below and cooled from above [TH5]. A well 
established feature of the dynamics of the system is the 
large-scale circulation (LSC). It plays an important role 
in natural phenomena, including convection in the Arctic 
ocean, in Earth's outer core, in the interior of gaseous gi- 
ant planets, and in the outer layer of the Sun. The prop- 
erties of the LSC have recently been intensively studied 
in experiments [4T-E4"] , numerical simulations [HO , 
and models [MED [H CIQ EH HHH3J. For a complete 
overview of all literature in which certain aspect of the 
LSC are studied we refer to the recent review of Aiders, 
Grossmann, and Lohse [T]. 

In experiments the LSC is measured by using thermis- 
tors that are embedded in the sidewall, see e.g. [33], or 
by using small thermistors that are placed in the flow 
at different azimuthal positions and different heights, see 
e.g. [15] . Since the LSC transports warm (cold) fluid 
from the bottom (top) plate up (down) the side wall, 
the thermistors can detect the location of the up-flow 
(down-flow) by showing a relatively high (low) tempera- 
ture. In addition, there have also been a number of direct 
measurements of the LSC by particle image velocimctry 
(PIV) and laser Doppler velocimetry (LDV) (see, for ex- 
ample Refs. [5J [7J [31]) that complement the thermistor 



method. 

In this paper we investigate the properties of the LSC 
with results from direct numerical simulations (DNS). 
Though DNSs are limited both in Ra and in duration, 
i.e. in number of LSC turnover times, to smaller values 
than the experimental analogues the advantage is that in 
contrast to the experiments the full spatial information 
is available. We will take advantage of this in order to 
verify the algorithms employed in experiments to iden- 
tify the LSC orientation based on a limited number of 
probes. We will introduce two existing methods to de- 
termine the LSC orientation over time. The first method 
[S] determines the LSC orientation from a cosine fit to the 
azimuthal temperature profile measured with the probes. 
This method has been extremely successful in revealing 
important features, i.e. azimuthal meandering, reversals, 
and cessations, of the LSC EJ HH H7J- In the second 
method the LSC orientation is determined by using a sec- 
ond order polynomial fit around the maximum and the 
minimum. Using this method the sloshing mode of the 
LSC, which is caused by an off-center motion of the LSC, 
was discovered [T5H2T] . We note that this can not be ob- 
tained by the cosine fit method as this method assumes 
that there is no off-center motion of the LSC. 

So far these procedures have been applied to exper- 
imental data where the data of 8 azimuthally equally 
spaced probes were available. The benefit of the simu- 
lations over the experiments is that it is easy to place 
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a large number of probes in the flow. Here we placed 
up to 64 azimuthally equally spaced numerical probes at 
different heights into the numerical RB sample. With an 
arrangement of up to 64 probes we can determine how the 
extracted information on the LSC and plume dynamics 
depends on the number of probes that is used. We also 
visualized the flow by movies that show the flow dynam- 
ics in horizontal or vertical planes, see the accompanying 
material [55] , 

The paper is organized as follows. In section [IT] we 
start with a discussion of the numerical method that has 
been used. This is followed in section [ill] by a discussion 
of the traditional methods that are used in experiments 
to determine the LSC orientation. Based on the results 
obtained by these methods we discuss the characteristics 
of the LSC in a T = 1 (section [TV} and T = 1/2 (sec- 
tion |V| sample. An important question for both high 
Ra number 36-38 and rotating RB convection [221 H3] 
is whether or not there is a single LSC present. There- 
fore we will discuss a new method in section \Vl\ which 
is based on the energy in the different Fourier modes of 
the azimuthal temperature and azimuthal vertical veloc- 
ity profile, to determine whether a single LSC is present. 
In section [VII| we will look into time resolved properties 
of the LSC. The general conclusions of this paper will be 



presented in section VIII 



II. NUMERICAL METHOD AND PROCEDURE 

The flow is solved by numerically integrating the three- 
dimensional unsteady Navier-Stokes equations within the 
Boussinesq approximation. For a detailed discussion of 
the numerical code we refer to Refs. [25l [33 [3S ED, • The 
flow is simulated in a cylindrical sample in order to keep 
the geometry identical to the one used in most experi- 
ments. In (551 [231 [371 13H] we have shown that simula- 
tion results obtained with this code agree excellent with 
experimental results. The cases we studied are based 
on the most common experimental setups that are avail- 
able, namely Pr 6.4 and Pr ss 0.7 at an aspect ratio 
r = D/L of 1 and 1/2. A detailed overview of the simu- 
lations can be found in table Q] 

In the simulations we placed up to 64 azimuthally 
equally spaced numerical probes per horizontal plane 
that provide simultaneous point-wise measurements of 
the temperature and the three velocity and vorticity com- 
ponents at the heights 0.25L, 0.50L, and 0.75L and a dis- 
tance 0.45D from the cylinder axis. Grid refinement tests 
were performed because the region close to the sidewall, 
where the LSC properties are sampled, is most sensitive 
from a resolution point of view |37j . The simulations are 
fully resolved according to the criteria of Stevens et al. 
[37] and the LSC properties we find do not depend on 
the grid resolution. Note that the azimuthal and radial 
number of grid points required to get the same resolution 
with respect to the turbulent length scales is less for the 
r = 1/2 cases than for the T = 1 cases because for the 



TABLE I. The columns from left to right indicate the follow- 
ing: Ra, Pr, F, the number of grid points in the azimuthal, 
radial and axial directions (Ng x 7V r x 7V Z ). The last two 
columns indicate the rel ativ e LSC strength determined using 
equation see section VI using the data of only 8 (5 m (8)) 
and all 64 (5 m (64)) probes in the horizontal midplane, re- 
spectively. 



Ra 


Pr 


r 


Ng x N r x N z 


5m (8) 


5 m (64) 


1 x 10 s 


6.4 


1.0 


257 x 129 x 257 


0.74 


0.70 


1 x 10 8 


6.4 


1.0 


385 x 193 x 385 


0.65 


0.68 


5 x 10 8 


6.4 


1.0 


257 x 129 x 257 


0.68 


0.73 


5 x 10 8 


6.4 


1.0 


385 x 193 x 385 


0.70 


0.76 


1 x 10 8 


0.7 


1/2 


193 x 65 x 257 


0.45 


0.57 


1 x 10 8 


0.7 


1/2 


257 x 97 x 385 


0.42 


0.54 


1 x 10 8 


6.4 


1/2 


193 x 65 x 257 


0.06 


0.27 



latter a larger horizontal extent has to be simulated. 



III. METHODS TO DETERMINE THE LSC 
ORIENTATION 

Following Ahlers and coworkers [TJ O H2], the orien- 
tation and strength of the LSC can be determined by 
fitting the function 



cos 



(1) 



to the temperatures recorded by the numerical probes at 
the height z = 0.5L and the azimuthal positions cf>i = 
2iir/n p , where n p indicates the number of probes. The 
three parameters 9 m , S m , and <p m are obtained from least 
square fits. Here S m is a measure of the temperature 
amplitude of the LSC and </> m is the azimuthal orientation 
of the LSC at midheight. The azimuthal average of the 
temperature at the horizontal mid-plane is given by 9 m . 
We calculated temperatures 9 t and 9b, orientations 4>t 
and 4>b, and amplitudes 5 t and St, for the top and bottom 
levels at z = 0.75L and z — 0.25L separately by the same 
method. For experiments [2 [HI [12] exactly this method 
is applied to the data of 8 equally spaced thermistors, see 
figure [T] 

A second method that is applied to determine the 
LSC orientation is to make a polynomial fit around 
the sensor that records the highest (lowest) temperature 
[191 120] . Finally, since we have a large number of numer- 
ical probes, i.e. up to 64 azimuthally equally spaced per 
horizontal level, we also determine the LSC orientation 
by finding the probe that records the highest (lowest) 
temperature. Again the LSC orientation based on these 
two methods is determined separately for all three levels, 
see figure [l] 

To distinguish the different methods we introduce a no- 



tation with a second index, i.e. 



^*,[c,p,m] 



. Here c indicates 



that 4> is determined using a cosine fit to the azimuthal 
temperature profile, p that <fr is determined using a poly- 
nomial fit around the maximum (or minimum) observed 
in the azimuthal profile, and m indicates that the maxi- 
mum (or minimum) observed in the azimuthal profile is 
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FIG. 1. Azimuthal temperature profile averaged over 4t/ at 
midheight for Ra = 5 x 10 8 , Pr = 6.4, and T = 1. The 
measured data from the numerical probes are indicated by 
the dots. The solid black line shows the cosine fit (eq. (TTj) ) 
to the data with offset m (dashed line) and amplitude dm- 
The black cross indicates the position of 4>m,c- The red cross 
indicates the position of the probe that records the highest 
temperature and thus the position of </>m >m . The blue solid 
line gives the polynomial fit around the maximum. Its max- 
imum, indicated by the blue cross, indicates the position of 
( t>m,p- The 4>{x) range that is considered in the polynomial 
fit is |7r. Panel a and b shows the procedure applied to the 
data of 8 and 64 equally spaced probes at midheight, respec- 
tively. The corresponding azimuthal vertical velocity profiles 
are given in figure [8] 



used. Because the numerical probes record both the tem- 
perature and vertical velocity component we use both to 
determine the LSC orientation. In the rest of the paper 
we will use the notation ^„ to indicate that the orien- 
tation is based on the azimuthal temperature profile and 
<fif* when it is based on the azimuthal vertical velocity 
profile. Hence <j) has three indices, i.e. <j)\ 2 , see tab' 
for a detailed overview. In the rest of the paper we wT 
indicate the numerical value of an index when this index 
is varied. 



TABLE II. The LSC orientation in this paper is determined 
from the azimuthal temperature or azimuthal vertical velocity 
profile with different methods. In addition, the LSC orienta- 
tion is also determined at different heights. As explained in 
the text this information is indicated in the notation <f>\ 2 . In 
this table the meaning of the symbols at the index locations 
1, 2, and 3, is summarized. 



Index 


Symbol 


Meaning 


1 


b 


Height 0.25L 


1 


m 


Height 0.50L 


1 


t 


Height 0.75L 


2 


c 


Cosine fit 


2 


V 


Polynomial fit 


2 


m 


Maximum 


3 


T 


Azimuthal temperature profile 


3 


w 


Azimuthal vertical velocity profile 



Because the numerically obtained azimuthal profiles 
can be very noisy, we applied a moving averaging filter 
to the data obtained from the numerical probes to elim- 
inate the detection of very small plumes. In experiments 
these events are not detected anyhow by thermistors, 
since these need time to react to temperature changes 
in the flow. We decided to apply a moving averaging fil- 
ter of 4tj, where 77 is defined with respect to the free-fall 
velocity Uf as 77 = L/Uf (L is the height of the sam- 
ple). Recently, Bailon-Cuba et al. [JT] showed the char- 
acteristic convective velocity U c of the LSC is approxi- 
mately Uf/5 for the parameter ranges Ra = 10 7 — 10 9 , 
T = 1/2 - 12, and Pr = 0.7. This means that the LSC 
turnover time t^sc defined as Tlsc — 2L/U C is about 10 



IV. RESULTS FOR T = 1 

Figure [2] shows the typical behavior of the LSC orien- 
tation for Pr = 6.4 in a T = 1 sample. The figure shows 
</>^ 2 an d 2- -"- n order to show that the LSC orientation 
based on the temperature and the vertical velocity data 
is almost identical we consider first a long averaging time 
(50r/). The results are displayed in figure^ and b. The 
similarity is expected since the LSC carries warm fluid 
from the bottom plate up the sidewall and vice versa. To 
show the influence of this long time averaging we show 
</>™ 2 based on the azimuthal vertical velocity profile af- 
ter averaging over 4t/, see figure [2ji and [2ji. When no 
time averaging is applied the graph looks similar, with 
some additional peaks due to very small plumes, to the 
one where the data is averaged over 4r/. Figure [3] shows 
that the LSC orientation can be determined more pre- 
cisely from the vertical velocity than from the tempera- 
ture, even for this relatively high Pr (i.e. small thermal 
diffusivity) . The difference between the LSC orientation 
determined from the vertical velocity and the tempera- 
ture seems to depend on several parameters such as the 
Ra and the Pr number, and the aspect ratio. However, 
at the moment we do not have enough numerical data 
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FIG. 2. The LSC orientation for Ra = 1 x 10 s , Pr = 6.4, 
and r = 1 on the 257 x 129 x 257 grid based on the data 
of 64 azimuthally equally spaced probes per horizontal level. 
The black, red, and blue line indicate 4>m,c 4>'rn,mi a ~d 4>m, P , 
respectively, a) <f>m,2 based on the azimuthal temperature 
profile averaged over 50r/ (w 5 LSC turnover times), b) 4>m2 
based on the azimuthal vertical velocity profile averaged over 
50r/, c) (j>m,2 based on the azimuthal vertical velocity profile 
averaged over 4t/. d) Enlargement of the graph shown in 
panel c to reveal more details. 



available to systematically study this. For T = 1 the ob- 
tained LSC orientation at midheight is the same when 
the data of only 8 or all 64 probes is used, but not when 
the data of only 2 probes is used. 

From figure [2] it becomes clear that the LSC orienta- 
tion, obtained by the polynomial and cosine fit, can differ 
significantly. This is due to the off-center motion of the 
LSC [TSHZE]. Figure [| shows = <f> max - 4> min , i.e. 
the difference between the orientation of the strongest 
up {(f>max) and down going motion ((j>min) obtained us- 
ing the polynomial fit, fluctuates around A<p = tt. To 
quantify the strength of these fluctuations we calculate 
((A^-tt) 2 ) 1 / 2 . The values based on the data of 64 (8) 

S robes per horizontal level for the case presented in figure 
are ((A^-Tr) 2 ) 1 / 2 ~ i. 09 (1.04), ((A^ p - tt) 2 ) 1 ^ « 
0.84(0.77), and ((A^ p - tt) 2 ) 1 / 2 w 1.11(1.01). Because 
the fluctuation strength is the same when the data of 8 
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FIG. 3. </>m c (solid black line) and <j> miC (red dashed line) for 
Ra = lx 10' 8 , Pr = 6.4, and r = 1 on the 257 x 129 x 257 grid, 
and based on the azimuthal vertical velocity and azimuthal 
temperature profile averaged over 4r/. Here the data of all 
64 azimuthally equally spaced probes at midheight have been 
used. 



and 64 equally spaced probes per horizontal level are con- 
sidered we conclude that the use of 8 probes is sufficient 
to capture these statistics in a T = 1 sample. In figure 
[5] it is shown that the LSC orientation, obtained using 
the polynomial fit method, can be different at the dif- 
ferent heights. We note that the same is observed when 
the LSC orientation based on a cosine fit is considered. 
This indicates that the LSC is not always flowing straight 
up and down, but is also moving in the azimuthal direc- 
tion, and thus performs twisting motions 6J. We note 
that some phenomena, like the drift of the LSC due to 
the Coriolis force, see Ref. [10J, can be analyzed using 
the LSC orientation based on the polynomial and cosine 
method, because in this case only the long term drift of 
the LSC is important. 

More remarkably is the presence of a specific frequency 
of approximately 25tj in the signal Cc~C.pi see figure 
[7j This specific frequency seems to be related to the fre- 
quency in which plumes are passing the horizontal mid- 
plain 42J, which in our case shows up in the frequency 
in which plumes are detected on the left and right hand 
side of the LSC orientation based on the cosine fit, see 
figure [6j This frequency might also be related to the low 
frequency mode found already in Ref. [43] or to the off- 
center oscillation of the LSC [55]. We find this frequency 
in our simulations for Pr = 6.4 in a T = 1 sample. We 
note that we also find this frequency when we do not ap- 
ply any time averaging on the data before we determine 
the LSC orientation. Actually one can already see this 
phenomenon in figure [2ji where 4>m,p fluctuates around 
Cc with a typical period of about 25r/. Further con- 
firmation is obtained when the temporal behavior of the 
data obtained from the 64 probes at midheight is shown 
in a movie, see the supplementary material 35 . The 
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FIG. 4. The red, black, and blue lines indicate A<f>f iP , A<f>^ p , 
and A0™ p , respectively for Ra = 1 x 10 8 , Pr = 6.4, and 
r = 1 on the 257 x 129 x 257 grid. Panel a and b show 
A</>™ p when the azimuthal vertical velocity profile is averaged 
over 50t/ and 4r/, respectively. The average of A0™ p is 7r 
and this is in agreement with a cosine fit. The deviations 
from this value are due to plumes and the off-center motion 
of the LSC. Here the data of all 64 azimuthally equally spaced 
probes per horizontal level have been used. 
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FIG. 5. The red, black, and blue lines indicate (j)f p , 4>m, P , an d 
4>™ p , respectively, for Ra = 1 x 10 s , Pr = 6.4, and T = 1 on 
the 257 x 129 x 257 grid. Panel a and b show <J>i iP based on 
the azimuthal profiles of the vertical velocity averaged over 
5O17 and 4r/, respectively. The difference between the LSC 
orientation at the different levels is due to plumes and the tor- 
sional motion of the LSC. Here the data of all 64 azimuthally 
equally spaced probes per horizontal level have been used. 



movie is based on the data obtained from the simulation 
on the 385 x 193 x 385 grid at Ra = 5 x 10 s with Pr = 6.4 
in a r = 1 sample. The movie shows a passing plume at 
t w 2412, see the snapshot of the movie in figure [HJd, 
which is indicated by the two peak structure around the 
maximum position obtained by the cosine fit. When this 




FIG. 6. Sketch of the top view on the RB sample. The dotted 
line indicates the orientation of the LSC based on the cosine 
fit and the red (up-pointing) and blue (down-pointing) arrows 
indicated the flow direction close to the bottom and top plate, 
respectively. Strong plumes are detected alternatingly on the 
left and right side of the LSC, see crosses, with a specific 
frequency, see also figure [2] 

plume is passed the double peak structure disappears. 
Subsequently, the LSC orientation based on the polyno- 
mial fit method is smoothly passing the LSC orientation 
based on the cosine fit. Thus showing an off-center mo- 
tion of the LSC. We note that having a large number 
of probes is essential to confirm these plume events. To 
show this we made a movie of the same time interval, 
see the supplementary material 35j and the snapshot in 
figure |8^, with the data of only 8 probes. From these 
movies and figure [8^ we conclude that it is impossible 
to see the double peak structure when the data of only 
8 probes is used and therefore we cannot distinguish the 
off-center motion of the LSC from passing plumes, when 
the plumes stretches over two rows of thermistors. 

In summary, in a T = 1 sample the obtained LSC ori- 
entation using a cosine and polynomial fit is the same 
when the data of only 8 or all 64 azimuthally equally 
spaced probes per horizontal level are used. Because the 
azimuthal temperature profile is only an indirect mea- 
sure of the LSC orientation we find that the LSC orien- 
tation can be determined with more precision from the 
azimuthal vertical velocity profile. In addition, we find 
that many of its interesting phenomena like azimuthal 
meandering, reversals, cessations, and the sloshing mode 
of the LSC look the same when the data of only 8 or 
all 64 azimuthally equally spaced probes per horizontal 
level are considered. However, only with the data of 64 
equally spaced probes per horizontal level the effect of 
passing plumes on the azimuthal temperature and az- 
imuthal vertical velocity profiles can be identified. 



V. RESULTS FOR T = 1/2 

We now come to the r = 1/2 case. For this geom- 
etry Verzicco and Camussi |25j have used the data ac- 
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FIG. 7. The power spectra of c 



shows a maximum 



around 25r/, which is indicated by the vertical dashed line. 
This frequency is caused by plumes that pass the LSC orienta- 
tion based on the cosine fit on the left and right side . The dif- 
ferent spectra are obtained from simulations for Pr = 6.4 in a 
r = 1 sample for different Ra and resolution, i.e. Ra = 1 X 10 s 
on a 257 x 129 x 257 (red) and on a 385 x 193 x 385 grid 
(black), and Ra = 5 x 10 s on a 257 x 129 x 257 (green) and 
a 385 x 193 x 385 (blue) grid. 



cessibility provided by DNS to show that the large scale 
circulation can be either in a single-roll state (SRS) or 
in a double-roll state (DRS). Subsequently, a model in- 
troduced by Stringano and Verzicco [35] found that the 
switching between the DRS and the SRS can influence 
the Nusselt number, and can thus be a possible reason 
for the bimodal behavior of Nusselt found in some exper- 
iments, see e.g. [33J |35] . 

The first extensive experimental study on the transi- 
tion between the SRS and the DRS was performed by 
Xi and Xia [16] . They studied flow mode transition in 
samples of aspect ratio 1, 1/2 and 1/3. Figure 10 of their 
paper [T5J shows the percentages of the time the flow 
spends in cither the SRS or the DRS. In addition, they 
showed that the heat transfer in the SRS is higher than 
in the DRS. Recently, Weiss & Ahlers [35J have also ex- 
perimentally investigated the switching between the SRS 
and the DRS, but now over a much larger range of Ra 
than Xi and Xia [16] . The conclusion of this work is sum- 
marized in figure 11 of Ref. [35] . That figure shows that 
the SRS is more dominant at higher Ra and smaller Pr. 
They note that they find a good agreement with the re- 
sults of Xi and Xia in the Ra number range in which both 
experiments overlap. In addition, the work of Xi and Xia 
[15] and Weiss and Ahlers [46] showed that the flow state 
can not always be defined as SRS or DRS. This is partic- 
ularly important for Ra <3x 10 9 . E.g., for Ra = 1 x 10 8 
and Pr — 4.38 the flow state is undefined for about 50% 
of the time according to Weiss and Ahlers [46] . 

In this paper we will qualitatively confirm these ex- 
perimental and numerical findings. We will restrict our- 
selves to this relatively low Ra number regime, namely 
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FIG. 8. Snapshot of the azimuthal vertical velocity profile 
averaged over 4r/ for Pr — 6.4 in a F = 1 sample at Ra = 
5 x 10 s . The symbols have the same meaning as indicated in 
figure [T] Panel a indicates the obtained profile using the data 
of 8 probes and panel b the profile as it is obtained from the 
data of 64 probes. The accompanying movies can be found in 
the supplementary material |35| . 



to Ra = 1 x 10 s : First, in order to directly compare with 
the T = 1 case of the previous section, second, because 
current high Ra number simulations |37[ 147] are too lim- 
ited in the number of turnover times to obtain sufficient 
statistics on the different flow states. 

Our simulations confirm Weiss & Ahlers's [35] finding 
that there is more disorder in a V = 1/2 sample than in 
a r = 1 sample. We find that in this low Ra number 
regime it is important to have a large number of probes 
to accurately represent the azimuthal temperature and 
azimuthal vertical velocity profiles to analyze the flow. 
In this regime we see that deflection of plumes due to 
corner rolls and other plumes happens more often than 
in a T = 1 sample. These effects can be observed in the 
supplementary movies for T = 1/2 at Ra = 1 x 10 8 and 
Pr = 0.7 and Pr — 6.4 [35]. In the movies that show 
horizontal cuts at the heights 0.25L, 0.50L, and 0.75L 



FIG. 9. The PDF of |<^ c - <p% c \ for Ra = 
for Pr = 6.4 and T = 1 is given in blue, 
line give the PDF for Pr = 6.4 and Pr - 
geometry, respectively. 



= 1 x 10 8 . The PDF 
The red and black 
: 0.7 in the T = 1/2 



one can see that the hot and cold regions interchange 
more often in a T = 1/2 sample than in a T = 1 sample. 

In order to determine whether there is a SRS or a DRS 
in this geometry we determined the probability density 



function (PDF) of the quantity |0£ 



IJ t,c\- 



When this 



PDF peaks around zero it means that the fluid is flow- 
ing straight up, and this indicates the presence of a SRS. 
When there is a DRS the PDF should peak around tt. 
Figure [9] shows that for Pr — 6.4 in a T — 1 sample the 
PDF peaks around zero and this confirms the dominance 
of the SRS in this geometry. For Pr = 0.7 in a T = 1/2 
sample the figure shows that there is a small peak around 
zero, but it is much less pronounced than for the T = 1 
case. For Pr = 6.4 in a V — 1/2 sample the figure shows 
that the PDF is nearly uniform, which would suggest that 
neither the SRS or the DRS is dominant. This agrees 
with the results of Weiss & Ahlers [46], since they find 
that for Pr = 4.38 and Ra < 3 X 10 9 the flow state is 



poorly defined and no state dominates. Figure 10 reveals 
the origin of these nearly uniform PDFs in the r = 1/2 
sample by showing three-dimensional visualizations of 
the regions where the vertical velocities are strongest. 
An analysis of these (and similar figures) revealed that 
the uniform PDF is due to the possibility of different flow 
states, i.e. the SRS (figure [Ml) with - <Pf c \ « 0, 
the DRS (figure [lOfc) with |0^ c - 4>™ c \ « 'tt, strong cor- 
ner flows (figure |10t>) which can have any |0™ c — <j)f c \, 
and strong plume and/or torsional motions (figure 10 i) 
which also can have any |</>™ c — (j>f c \. 

To better compare our results with the results of Weiss 
& Ahlers [33] we also determined the PDF of the quan- 
tity <\>™i t c — 4>m,ci see figure |TTj Here </>^, t c means that 
both (j)™ c and (j>f c are compared with 0™ c to construct 
the PDF. A similar PDF for higher Ra, based on experi- 
mental data, can be found in figure 21 of Weiss & Ahlers 
[4"S] . Figure 11 shows this PDF based on the data of 





FIG. 10. Snapshot of the regions with the strongest, i.e. 
w ^ 0.35wmax, up (red) and down-flow (blue) for Ra = lx 10 s 
and T = 1/2. a) Pr = 0.7, SRS, b) Pr = 0.7, very 
strong corner flows, c) Pr = 6.4, DRS, d) Pr — 6.4, strong 
plume/torsional motion. Supplementary movies that offer a 
full three dimensional view are provided in the supplementary 
material [35] . 



only 8 and all 64 equally spaced probes per horizontal 
level for the cases of figure [9] First of all, figure [TT] shows 
that the DRS is much more pronounced for Pr = 6.4 
than for Pr — 0.7, which is in agreement with the re- 
sults of Weiss and Ahlers [46] and the results shown in 
figure [9] Second, a comparison of the PDFs in figure [TT] 
reveals that the PDFs based on the data of 8 and all 64 
equally spaced probes are almost identical. The largest 
difference is visible for Pr = 6.4 in a T = 1 sample. Here 
the PDF based on the data of 64 equally spaced probes 
is more confined in the SRS region than the one based 
on 8 equally spaced probes, because the larger number 
of probes reduces 'random' fluctuations in the obtained 
LSC orientation. 

Furthermore, and this is the main issue here, for this 
relatively low Ra sometimes the data of three rows of 8 
thermistors can give misleading hints for the actual flow 




A(C/ t -4> )(rad) 

VT b/t,c T m,c' v ' 



FIG. 11. The PDF of <ffi /%c - 4C,c for Ra = 1 x 10 s , see 
text for details. The dashed and solid lines indicate the PDF 
based on the data of only 8 or all 64 probes equally spaced 
probes per horizontal level, respectively. The vertical dashed 
lines indicate the definition of the SRS according to Xi and 
Xia US and Weiss & Ahlers ED. 
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FIG. 12. Red (top), black (middle), and blue (bottom) data 
indicate the azimuthal temperature profile at the heights 
0.25L, 0.50L, and 0.75L, respectively, for the flow state in- 
dicated in figure |10b . The dots indicate the data of 8 az- 
imuthally equally spaced probes per horizontal plane. Based 
on the cosine fits (solid lines) obtained from this probe data 
this state would be identified as a DRS, while figure |Tup shows 
that it is in fact a SRS with a strong corner flow. The dashed 
lines show that a higher azimuthal resolution can resolve the 
double peak structure of the corner flow and the main flow. 
The temperatures at 0.25L (0.75L) have been raised (lowered) 
by 0.1 for clarity. 



state of the system due to the effect of passing plumes and 
corner rolls on the azimuthal profiles. The importance of 
these effects is shown in figure [12] This figure shows that 
based on the data of 8 azimuthally equally spaced probes 
at each horizontal level the flow structure shown in figure 
|10b, which is a SRS with a strong corner flow, would be 



identified as a DRS. In addition, the figure shows that the 
full azimuthal temperature profiles reveal that a higher 
azimuthal resolution (e.g., 64 instead of only 8 probes) 
could resolve the double peak structure (originating from 
the coexistence of the SRS with some corner flow), which 
makes it possible to distinguish this state from a real 
DRS. 

We note that the experiments of Weiss & Ahlers [4"6] 
show that for Pr = 4.38 and Ra < 3 x 10 9 the flow 
state is undefined for about 50% of the time. In the ex- 
periments the SRS only establishes itself for higher Ra, 
where the flow state is much better defined, as the area 
of the corner flow becomes smaller and smaller [23]. It 
is therefore likely that the examples presented here are 
primarily important in the low Ra number regime in- 
vestigated here and are much less common in the higher 
Ra number regime considered in most experiments. We 
stress however that one does not know a priori whether 
one is in a difficult case or in a case where the method 
works fine. Therefore we feel that it is important to re- 
alize that there are some limitations in the methods to 
determine the flow state from three rows of 8 equally 
spaced thermistors. An example of such a difficult case 
is the breakdown of the LSC when a strong rotation is 
applied around the cylinder axis [221 [23] . 



VI. THE RELATIVE LSC STRENGTH 

The instantaneous temperature profiles in figure |12| 
show that the LSC is not always clearly present, meaning 
that the cosine fit becomes poor. In order to study this 
in more detail, we introduce the concept of the relative 
LSC strength, which we define as the strength of the LSC, 
i.e. the cosine mode, with respect to the strength of the 
plumes and turbulent fluctuations, i.e. the fluctuations 
around the cosine fit. 

A convenient way to define a relative LSC strength that 
is independent of the number of probes is to determine 
the energy in the different Fourier modes of the azimuthal 
temperature or the azimuthal vertical velocity profile. In 
order to obtain a relative LSC strength Sk with a number 
between and 1 we normalize the energy in the first 
Fourier mode, i.e. the cosine mode, as follows: 



S k = Max 




1 

N 



.0 



(2) 



Here ^ t e E\ indicates the sum of the energy in the first 
Fourier mode over time, i.e. from the beginning of the 
simulation t = % to the end of the simulation t = t e , 
Y] t e E to t the sum of the total energy in all Fourier modes 
over time, and N the total number of Fourier modes that 
can be determined. Note that the factor 1/N is only 
used to subtract the part of energy coming from the equi- 
distribution among all the modes. The number of Fourier 
modes A" that can be determined is n p /2, where n p — 2 Z 
is the number of probes and i is an integer. The subscript 
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k (= b, m, t) indicates the height level in the sample at 
which the relative LSC strength is determined. 

Here we will only determine Sk for the middle row 
of thermistors, thus only S m . We note that the LSC 
amplitude 6 m and the phase <j> m of the first Fourier mode 
are identical to the ones determined when using the least 
square fit in equation ([!]). We moreover note that the 
quantity Sk is not time resolved and thus does not allow 
to study the switching between the SRS and DRS [4*o] . 
Results on time-resolved relative LSC strengths will be 
presented in the next section. 

The normalization in eq. ^ has been chosen such that 
the relative LSC strength S m always has a value between 
and 1: The value 1 indicates that the azimuthal profile 
is a pure cosine profile, which is a signature of the LSC 
according to [5] , and the value indicates that the magni- 
tude of the cosine mode is equal to (or weaker than) the 
value expected from a random noise signal. S m ^> 0.5 
indicates that the cosine fit on average is a reasonable 
approximation of the data, as then most energy in the 
signal is in the first Fourier mode, i.e. the cosine mode. 
In contrast, S m <C 0.5 indicates that most energy is in the 
higher Fourier modes, meaning that the application of a 
cosine fit to the data becomes questionable. Somehow 
arbitrary, we define the SRS as dominant once S m > 0.5. 

To further demonstrate these properties of Sk, we an- 
alyze the test function 



0i(j) = A cos cos(& + j/100) + f(j). 



(3) 



Here <pi = 2iir/n p and f(j) is a set of Gaussian dis- 
tributed random numbers with a certain standard de- 
viation (J 2 ) 1 / 2 , the "time" j is varied between 1 and 
125000, and the noise is supposed to model the effect of 
turbulent fluctuations around the cosine fit. The result 
for the relative strength Sk of the LSC as function of the 
noise level is shown in figure 13 Indeed, 1 > Sk > and 
Sk > 0.5 only if the noise is smaller than the amplitude, 
(f 2 ) 1/2 /A cos «l. 

We now come back to the numerical data of our sim- 
ulations. For each simulation we calculated the relative 
LSC strength S m based on the instantaneous azimuthal 
vertical velocity profiles at midheight based on the data 
of either only 8 or all 64 probes at the horizontal mid- 
plane. We show the results for two aspect ratios, two 
Rayleigh numbers, and two Prandtl numbers in table [T] 
For r = 1 the flow is clearly in the SRS, as indicated by 
S m > 0.65. For r = 1/2 the SRS is less dominant, in 
particular for the large Prandtl number case Pr = 6.4, 
where the SRS hardly occurs. This confirms our conclu- 
sions of the previous section where we stated that the 
SRS is not always present in the T = 1/2 sample. 

Another method used in the literature 9, 46] to iden- 
tify the flow pattern is to average the azimuthal profile 
with respect to the LSC orientation. This method as- 
sumes that the deviations from the cosine fit are due to 
turbulent fluctuations. It will work for cases for which 
S m is large. However, when S m is low, this method gives 




COS 



FIG. 13. The relative strength of the LSC Sk versus the rel- 
ative noise level (f 2 ) 1 ^ 2 /A coa , calculated from the stochastic 



model (equation 
is defined to brea 



3f ). For strong noise Sk < 0.5 and the SRS 
c down. 



< 
< 




misleading results as is shown in figure 14 where we apply 



FIG. 14. The azimuthal vertical velocity profile after aver- 
aging with respect to the orientation of the LSC, which was 
obtained using a cosine fit, is indicated by the dots. The 
method has been applied to the data of all 64 probes in the 
horizontal midplane obtained in the simulation at Pr — 6.4, 
Ra = 1 x 10 8 , and F = 1/2. The solid line indicates a pure 
cosine profile. 



it to above discussed data for Pr = 6.4 in the T = 1/2 
sample, where the relative LSC strength using the data 
of all 64 probes in the horizontal level is S m = 0.27 (see 
table []]): The average profile obtained with the method 
of refs. [9l [46] cannot be distinguished from the ones ob- 
tained from the cases where S m is larger and therefore 
the method cannot be used to claim the dominance of 
the SRS. 
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FIG. 15. PDFs of the relative LSC strength at midheight for 
Ra = 1x10 s , Pr = 6.4, and T — 0.5. The red (-§ < S m (i) < 
1) and blue (— < S m (t) < 1) data points indicate the PDF 
based on the azimuthal vertical velocity profile sampled with 
only 8 and all 64 probes, respectively. The colored dashed 
lines give the corresponding PDFs based on a purely random 
signal. The vertical dashed lines at S m = 0.5 indicates the 
region 0.5 < S m < 1.0 where the first Fourier modes contains 
at least 50% of the energy of the azimuthal vertical velocity 
profile. 



VII. TIME RESOLVED RELATIVE LSC 
STRENGTH 

The relative LSC strength can also be calculated on in- 
stantaneous azimuthal temperature and vertical velocity 
profiles. We demonstrate this by constructing the PDF 
of the relative LSC strength at midheight S m (t) for in- 
stantaneous azimuthal temperature and vertical velocity 
profiles using the definition 



Etot{t) 



1 

N 



I 1- 



(4) 



The only difference with respect to equation ^ is that 
Exit) and E tot (t) and thus S m {t) are time dependent 
and that we have dropped the criterion that the minimal 
value should be 0. The reason for the latter is that we 
want to prevent strange jumps in the PDFs of S m at 0. 
For 8 probes we have — A < S m {t) < 1 and for 64 probes 
we have — ^ < S m (t) < 1. Note that these intervals 
follow directly from equation [4] by filling in the limiting 
cases, i.e. Ex(t) / E tot (t) — and Ei(t) / E tot (t) = 1, and 
by using the relation N — n p /2, where n p is the number 
of probes, see previous section. 

Figure [15] shows the PDFs for S m based on the mea- 
surement of the vertical velocity by only 8 or all 64 
equally spaced probes at midheight in the simulation 
with Ra = 1 x 10 8 , Pr = 6.4, and in the T — 0.5 sample. 
In addition, for both cases we determined the distribu- 
tion of S m (t) for a random signal. When we compare the 
PDFs based on the simulation data with the distributions 



obtained for random signals, we see that the difference is 
relatively small when the data of only 8 probes is consid- 
ered, whereas a much larger difference is obtained when 
the data of all 64 probes is taken into account. Figure 
15 thus clearly shows that it can be beneficial to use 
more than 8 azimuthally equally spaced probes, since it 
becomes much easier to show that certain events are sta- 
tistically relevant. 



VIII. CONCLUSIONS 

We studied the LSC dynamics in DNS simulations by 
investigating the azimuthal temperature and vertical ve- 
locity profiles obtained from 64 equally spaced numerical 
probes at three different heights. For Pr — 6.4 in a 
r = 1 sample we find that the azimuthal profile is well 
presented with the data of 8 numerical probes, a number 
normally used in experiments, as the LSC orientation ob- 
tained by a cosine fit is the same when the data of 8 and 
64 probes is considered. We find that the improved az- 
imuthal resolution (64 instead of 8 probes in the sidewall 
at one height) can reveal the effect of the plumes. 

In agreement with the findings of Xi and Xia |16j and 
Weiss and Ahlers jUj we find that there is more disorder 
present in the V = 1/2 sample than in the T = 1 sam- 
ple. For the T = 1/2 case we also show that when the 
azimuthal temperature profile is only determined from 
8 probes per horizontal level a SRS can erroneously be 
identified as a DRS, because the azimuthal resolution is 
too small to distinguish between the structure of the cor- 
ner flow and the main roll. Here we again stress that 
this result is obtained in the relatively low Ra number 
regime. The experimental results of Weiss & Ahlers [4"S"] 
show that for the low Ra number regime the flow state 
is undefined for about 50% of the time. For T = 1/2 the 
SRS only establishes itself for higher Ra, where the flow 
state is much better defined [46]. It is therefore likely 
that the examples presented here are primarily impor- 
tant in this low Ra number regime investigated here and 
are much less common in the higher Ra number regime. 

We quantified the LSC strength relative to the turbu- 
lent fluctuations by determining the ratio between the 
energy in the first Fourier mode and the energy in all 
Fourier modes of the azimuthal temperature and az- 
imuthal vertical velocity profiles. We find that the rela- 
tive LSC strength at Ra = 1 x 10 s is considerably lower 
in the T = 1/2 sample than in the T = 1 sample, i.e. that 
the SRS is much less pronounced in the r = 1/2 sample 
than in the T = 1 sample. This determination of the 
relative LSC strength can be applied directly to avail- 
able experimental data to determine whether the SRS 
is present in high Ra number thermal convection and in 
rotating RB convection. 
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